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Abstract 

A new computational method for solving the nucleon-deuteron breakup scattering problem 
has been applied to study the elastic neutron- and proton-deuteron scattering on the basis of 
the configuration-space Faddeev-Noyes-Noble-Merkuriev equations. This method is based on the 
spline-decomposition in the angular variable and on a generalization of the Numerov method for the 
hyperradius. The Merkuriev-Gignoux-Laverne approach has been generalized for arbitrary nucleon- 
nucleon potentials and with an arbitrary number of partial waves. The nucleon-deuteron observ- 
ables at the incident nucleon energy 3 MeV have been calculated using the charge-independent 
AV14 nucleon-nucleon potential including the Coulomb force for the proton-deuteron scattering. 
Results have been compared with those of other authors and with experimental proton-deuteron 
scattering data. 

PACS numbers: 21.45. +v,11.80.Jy,25.45.De 
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I. INTRODUCTION 



There is an impressive amount of nucleon-deuteron scattering data: proton-deuteron 
and neutron-deuteron elastic and breakup data: total, partial and differential cross sections 
and spin observables involving nucleon and deuteron beams. The data are compared with 
the rigorous three-body theory: Faddeev-equations-based theory using as input realistic 
high-precision nucleon-nucleon potentials, and including model three- nucleon forces [l]. In 
some calculations Coulomb force has been included Nucleon-nucleon (NN) potentials 
used in rigorous three-nucleon (3N) calculations describe the NN database with x^/degree of 
freedom approximately equal to one. These are AV18 [3], CD-Bonn jj] and several Nijmegen 

6|. Among three nucleon forces (3NF) are Tucson- 



potentials j5| and to a lesser degree AV14 
Melbourne and its various modifications 



7|, and Urbana potentials jsj. Based on the chiral 



effective field theory (EFT) NN and 3N potential have been developed {q] and they have 
been used in a rigorous 3N calculations 10| ■ A local version of the effective field theory at 
next-to-next to leading order labeled N^LO is given in ref. 

In spite of this enormous progress in the three-nucleon studies, there are several important 
cases where the rigorous three-nucleon calculations have failed to explain the data [12'] and 
these discrepancies are established with very high precision. Among the most important 
discrepancies are the Ay puzzle in nucleon-deuteron (Nd) elastic scattering [l3|, the star 
configuration in the Nd breakup reaction [l^, quasi- free scattering (QFS) cross section |l| 
and the nd backward angle scattering at energies between 50 and 100 MeV ISj. Some three- 
nucleon data show clear evidence for the 3NF, but some are in better agreement with the 
calculation if the 3NF are not included. High precision realistic potentials (Nijmegen, Bonn, 
Paris, Urbana) are not phase equivalent and they predict different triton binding energies, 
they have different short range potentials and some differ conceptually. It is hoped that 
EFT will give an answer, but it is still unclear. 

There are more 3N data involving charged particles and therefore, calculations rigorously 
including electromagnetic interactions are of paramount importance. The pd scattering has 
been studied by using hyperspherical harmonic method and Kohn Variational Principle 1^ 
and by using the screening and renormalization procedure 17l|. At 3 MeV calculations have 
been done using high precision realistic potentials and 3NF [l8|, while at energies above the 
threshold for the deuteron breakup only calculations using screening and renormalization 



have been done. The screening method cannot be apphed to energies below 1 MeV and this 
is a serious hmitation. 

In this article we present the development of an alternative method for studying the 
proton-deuteron (pd) system based on the direct numerical solution of the Faddeev-Noyes- 
Noble-Merkuriev (FNNM) equations in configuration space. This approach was initiated by 
Merkuriev et al. (MGL) 19| who derived general formulae for nd breakup scattering. This 
method has been originally applied to study nd and pd elastic and breakup scattering but 
limited only to nuclear S-waves interaction and to simple NN potentials |20j. In the present 
work we generalize the MGL approach to any high precision realistic potential for both nd 
and pd for elastic processes. 

This paper is organized as follows: in section 2 we describe a calculation in configuration 
space starting with the general formalism in subsection 2.1, followed by Numerov method 
in subsection 2.2. Our novel method for solution is given in subsection 2.3. Our results 
are presented in section 3. Comparisons of our results with the previous calculations and 
with the data are discussed in section 4. Finally, our summary and conclusion are given in 
section 5. 



II. THREE-NUCLEON FADDEEV CALCULATION IN CONFIGURATION SPACE 
- OUR NEW COMPUTATIONAL METHOD 

A. Formalism 

The starting point for studying interactions between nucleons in three-body systems is 
the solution of the Schrodinger equation H'^ = E^^ for nuclear Hamiltonian such as 

^ = -|:Ev? + K + E^.^ (+ E ^.h), (1) 

i=l j<k j<k<l 

where Vc and Vjk are the Coulomb and nuclear potentials, respectively. In this study we 
neglected by three-nucleon forces Vjki- 
Writing the total wave function as 

^ = $i + $2 + $3 = (l + P+ + p-)$i, (2) 
the Schrodinger equation for three identical particles can be reduced into a single Faddeev 
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equation, which in Jacobi's vectors xi,yi has the form 



(3) 



where the operators are the cychc permutation operators for the three particles which 
interchange any pair of nucleons (P"*" : 123 — 231, : 123 — t- 321). The Coulomb potential 
has the following form: 



me 



a \^a\ 2 



n 



2 J 



(4) 



where e^=1.44 MeV -fm and h'^ /m=4:lA7 MeV-fm^. The sum runs over a =1,2,3 for the 
three possible pairs and the product of the isospin projection operators runs over the indices 
i of the particles belonging to the pair a. As independent coordinates, we take the Jacobi 
vectors Xa^y^. For the pair q;=1, they are related to particle coordinates by the formulas: 



xi = r2 - rs. 



(5) 



for a=2,3 one has to make cyclic permutations of the indexes in Eq.(5). The Jacobi vectors 
with different a's are linearly related by the orthogonal transformation 



^x 



CafS Safi 
-Sal3 Ca/3 



where 

CafS 



(r 

Xp 



^2 , c;-2 



(6) 



(M-m„)(M-m,)' ^ ^ sgn(/3-a)v'l-C^„ M = ^m^. (7) 

To perform numerical calculations for arbitrary nuclear potential, we use MGL approach 



19| . For pd scattering the FNNM equations for partials components can be written in the 



following form (here we omit the index 1): 



E + + dl) - y) $Ao.^o,M„(^^ ^ ^^^^^ + + v-c^\ 



.Mo-- 



(8) 



x$J°'^°'^"(x, y) + v.p{x) [^f^'''''\x, y) + S\ du g^iy/^. n)$^-^-^'^o(a;', y') . 

Here Greek subindexes denote state quantum numbers: a = {/, a,j, s, A, t, T}, where /, a, j 
and t are the orbital, spin, total angular momenta and isospin of a pair of nucleons, A is the 
orbital momentum of the third nucleon relative to the c.m.s. of a pair nucleons, and s is the 



total "spin" (s = l/2+j). M = A + sis the total three-particle angular momentum, and 
the value of total isospin is T. In Eqs.(l8]) vi and coefficients c'^^o^ depending on quantum 
state numbers of channel combined with are matrix elements of the Coulomb potential 
projected onto the MGL basis. For given a and /3 summation over r in Eqs.(IH]) is finite. 
If a = {lajsXtT} and /3 = {I'a'j's'X't'T'} then values of r are restricted by the following 
inequahty: 

max(|/ - /'I, A - A'l) < r < min(/ + A + A') 

. This means that for a chosen set of basic states, Eqs.(IH]) take into account the Coulomb 
interaction exactly (although the latter has been expanded in partial waves). 

The geometrical function gis^{x,y,u) is the representative of the permutation operator 
P+ + P- in MGL basis 3: 



(-1) 



ga'a(y/x, U) = ga'a{0, u) = Qa'aiO, 9') 

A+A'+j+j'[(2J + 1)(2J' + l)(2s + l)(2s' + 5](2S + 1)(2L + 1) 

LS 



I a J \\ I' a' J' \\ \ I L \\ X' V L 
1/2 s 5 J I 1/2 s' S \ \s M s \ \s M s' 

< xf/2a'Vl/2,t'\P^\X'!/2aVl/2,t > I' XliV / ^ , u) ■ (9) 

Function h is the representative of the permutation operator + P in the X + I = L 
basis: 

-^(-l)'+^i^^±iM±l)[(2A)!(2/)!(2A' + l)(2r + 1)]V^ 



x'y' 



E(-1)^2A;+1)P,H 

k=0 Xi+X2=X,li+l2=l 



yXi+h r^X2+h 



-r 



[(2Ai)!(2A2)!(2/i)!(2/2)!]V2 J^}^^" + ^^^^Z" + 1) 




(10) 



The index k runs from zero to (A' + /' + A + The (...) are the 3j symbols: 

•^1 -^2 J3 I _ / ^ yj^+rn3+2ji ^ fijam^ 

\ ^' /9„- I 1 ^ii-mij2-m2- 

mi 1712 ^3 / V^J3 -r -L 

The centrifugal potential is 

n^,i{i + i) A(A + i) 



m 



+ 



^2 ^2 



(11) 



and nucleon-nucleon potentials are Vaa'{.x) =< a\v{-x)\a' >= SxySss'Sf^a'^jj'vf/ , where v^/ 
are the potential representatives in the two-body basis yfj'' (x) (most often abbreviated as 

The set of partial differential equation Eqs.([8]) must be solved for functions satisfying the 
regularity conditions 

^AosoAfo^O, 6) = <l>^«"»*^o(p, 0) = <l>^«"o^^»(p, 7r/2) = (12) 

The asymptotic conditions for pd elastic scattering has the following form 2l[ |: 



(13) 

X finite, y ^ oo, 

where =argr(A + 1 + iu) is the Coulomb phase and u is equal n/{\/3q), tpi is I — th 
component of deuteron wave function (Z =0,2), and F'^ and are the regular and irregular 
Coulomb functions, respectively. 
The S-matrix is defined as follows 

S^slso = hxoSssMie''''^^ + 2ialXs,. (14) 
At energies below threshold the S'-matrix is unitary and may be presented as 

S = e^^^, 

where A is the Hermitean matrix of scattering phases. From ( !T3|) we then find that the 
matrix of partial elastic amplitudes a has the structure 



where A*^ is a diagonal matrix of Coulomb phases and A is the Hermitean matrix of scattering 
phases. The phases 6 = A—A^ are the contribution to the scattering phase due to the nuclear 
interaction. 

To simplify the numerical solution the FNNM equations, we write down Eqs.([H]) in the 
polar coordinate system (p^ = x"^ + y"^ and tan 6* = y/x): 



+ Y,v^p{p,9)[u^p^^-''\p,9) + f^duY.gp,{o,u,e\e,u))u^^^-^'^\p,e'] 



(16) 



/3 

where 

1 3 

cos^ 6 (u, 6) = - cos^ 6 cos ^ sin ^ ■ m + - sin^ 9, (17) 

V ' / ^ 2 4 

and the first derivative in the radius is eliminated by the substitution U = p~^/^$. In Eq. 
f fT6|) Qa/3 is the overall matrix sum of the Coulomb potential. 

In the case of neutron-deuteron elastic scattering one has to set the "charge" n equal to 
zero. This leads to equality to zero of the Coulomb phases A^, and the Coulomb functions 

and are reduced to the regularized spherical Bessel functions j\ and — ^a, respectively. 



B. Numerov method 

Modification of the Numerov method for the set of the differential equations ( 1161) does 
not present any difficulties in principle. As is well known, the Numerov method is an 
efficient algorithm for solving second-order differential equations. The important feature 
of the equations for the application of Numerov's method is that the first derivative has 
to be absent. The aim of this method is to improve the accuracy of the finite-difference 
approximation for the second derivative. Starting from the Taylor expansion truncated after 
the sixth derivative for two points adjacent to x„, that is for Xn-i and one sums these 
two expansions to give a new computational formula that includes the fourth derivative. 
This derivative can be found by straightforward differentiation of the second derivative from 



the initial second-order differential equation (see the details in 22(|). For brevity, we omit 



the corresponding derivation and present only the final formula of Numerov's method for 



the FNNM equations (omitting the upper indices AoSqMo): 

+ EK«(Pj,9)-Api4fl(Pj,9))(Cfl(Pj-i.e)+i:/_*<i9V,(9.9')f^7(ft-i,C')) 

[l^P) Pj Pj p Pj Pj 

r!, 12 2/\pT^{e). ' Ap 

\L\p) Pj Pj p Pj Pj 

+ j:Mpve) + Apv'^i,ip,,e))iu^ip,^,,e) + Y: ^ de'gp,{e,9')u,{p,^,,e')), 

n ^. J 6 



(18) 



/3 

where 



^ 9^^2 cos2 sin^^^ 4' 



In Eq. f lT8|) pj is the j — t/i current point for hyperradius p G (0, Rmax) in the radial grid 
(j = 1, 2, . . . , Np), Apj is the radial step-interval. 

To ensure the accuracy of order {A6)^ for the approximation in the angular variable, 



Hermitian splines of the fifth degree have been used (see Ref. 23|). These splines are local 
and each spline Scri{x) is defined for x belonging to two adjacent subintervals and 
[xjjXj+i]. Their analytical form is fixed by the following smoothness conditions: 

Sai{xi^i) = 0, Sai{xi+i) = 0, a = 0, 1, 2, (19) 

and 

Soiixi) = 1, 5'oi(xi) = 0, 5'oi(xi) = 0, 

Su{xi) = 0, S[,ix,) = 1, SUx^) = 0, (20) 

S2i{Xi) = 0, S'2i{Xi) = 0, S'^iiXi) = 1. 

Expansion of the Faddeev component into basis of the Hermitian splines has the following 
form: 

2 Ne+1 

u^{p,0) = y: E S.,i9)C:,ip)^ (21) 

(7=0 j=0 

where A'e + 1 is the number of internal subintervals for the angular variable 6 G [0, 7r/2]. 

To reduce the resulting equation (ITSl) to an algebraic problem, one should explicitly 
calculate the derivatives of NN potentials Va/3{p, 0) with respect to p and the second derivates 
of splines S„j{6) with respect to 6. It is convenient to express the second derivative of 
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component Ua with respect to 9 through Ua itself using Eg. ( 121]) . Upon substituting the sphne 
expansion fl2T]) and expression for its second derivative into Eqs. flT8|) . we use a collocation 
procedure with three Gaussian quadrature points per subinterval. As the number of internal 
breakpoints for angular variable 6 is equal to Ng, the basis of quintic splines consists of SNg+G 
functions. Three of them should be excluded using the last two regularity conditions from 
(fT2|) and continuity of the first derivative in 6 of the Faddeev component at either 6 = 
OT 6 = 7t/2, as the collocation procedure yields 3Ng + 3 equations. Finally Eqs. ffTSj) for the 
Faddeev components are to be written as the following matrix equation: 

A^U^ + G^U2 =0, 
BjU.^i +A,Uj + G,U,+,=0, j = 2,...N,-l, (22) 

^NpUwp-l + An^Unp = -GnpUnp+1- 

In this equation vector Uk = U{pk) has dimension A^j„ and matrices A, G have dimension 
Nin X Nin where Nin = x A^^, and is the number of partial waves and A^^^ = 3Ng + 3 
is the number of collocation points in the angular variable 6. 

C. The novel method of solution 

To derive equations for calculation of elastic Nd amplitudes, the method of partial inver- 
sion I22I has been applied. We write down Eq. (l22i) in a matrix form: 

{D*U), = -6,NpGNpUNp+i. (23) 

Here matrix D is of dimension NpNin x NpNin, and Np is the number of breakpoints in the 
hyperradius p. The form of this equation results from keeping the incoming wave in the 
asymptotic condition (|T3l) . As a consequence, the right hand part of Eq.f l23l) has a single 
nonzero term marked with index Ap + 1. Sparse (tri-block-diagonal) structure of matrix D 
optimizes considerably the inversion problem. 

Hyperradius pNp+i = Rmax, where Rmax is the cutoff radius at which the asymptotic 
condition Eq. fll3p are implemented. By formal inversion of the matrix D in Eq. fl23p . the 
solution of the problem may be written in the following form: 

= -DJ^GmpUnp+,, J = 1, 2....Np. (24) 
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In Eqs.f l24p one should consider the last component of vector U: 



- -D^^^^GnpUnp+i- (25) 

Provided Rmax is large enough, the vector U^p on the left side of Eq. (25) may be replaced 
by the corresponding vector obtained by evaluating Eq. f[T51) at the radius p = p^^. As a 
result in the case M > 3/2 we obtain three linear equations for the unknown amplitudes 



Mo 

3 



"AsAoso- 



E<°-v^ = F^ J = 1,2,3. (26) 

i=l 

For Mq= 1/2 the indices run over i,j=l,2. In these equations indices i,j number the 
asymptotic values of pairs (As), and vectors v, F are known quantities. For the sake of 
brevity, we do not display here the explicit form of them. As Rmax ~^ oo the set of equations 
( 126|) has a set of constants solution. At finite Rmax its solution is a vector a with 

generally different components corresponding to different angles. 

For each value of j linear equation ( 12 6 p is over determined, since the number of equations 
is Nin and the number of unknowns is 3. Therefore it is natural to use the least-squares 



method (LSM) as was proposed in 22|. According to LSM one has to minimize the following 
functional 

115^0^^° ■ V-Fi' = min. (27) 

i=l 

Differentiating this expression with respect to Rea^" and Im aj^° we obtain three (two) 
sets of liner complex equations of dimension 3 x 3 (2 x 2 for Mo=l/2), respectively. 



^a,f°-(v*,V) = (v*,F^), J = 1,2,3, fc = l,2,3, (28) 



3 

i=l 

where (^*, /) = Yjiii " fi is an ordinary scalar product. Now calculation of amplitudes 
is trivial task. 



D. Observables 

To calculate observables for elastic scattering of nucleon from deuteron in the direction 
q' (initial direction q is along the z-axis), one has to derive the equation for the elastic 
amplitude as a function of scattering angle. Omitting this derivation, we represent the final 
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expression for this amplitude in MGL basis: 



a'+J'^\Q,su:,+J-,^l/2a'lJL^l/2(7AjJ>^'M:,-a'^-Ji\^. J" 



-J's', 



M X's' As 



/2A 



47r 



(29) 



M 



with Mz = CTz + Jz- 

In Eq. fl29l) a'a'^{a, az) and J' J'^JJz) are spin and its projection for incoming (scattered) 
nucleon, and the deuteron in the rest (scattered deuteron), respectively. Thus, the nuclear 
part of the elastic amplitude is a (2 x 2) (g) (3 x 3) matrix in the spin states of nucleon and 
deuteron, depending on the spherical angles and 0. 

The situation is a little bit more complicated with the elastic scattering of proton from 
deuteron, since apart from the nuclear part the elastic amplitude also contains the pure 
Coulomb part. Thus in the matrix notation the resulting amplitude is to be sum of two 
amplitudes: 

a*"* = a + a\ (30) 

where a is the nuclear part of the same form as for the nd case and cf is the Coulomb part 
which is a unit matrix in spin states (this term does not change spins and depends only on 



«a^,J^,a,,J,(q) 



The amplitude is as follows 



n 



OTT qVSsm'^{e/2) 



(31) 



(32) 



The parameter v is defined by the ratio v = g is the wave vector of proton, the parameter 
n is given in Eq. and t]'^ = arg r(l + iu). 

The spin observable formulas can be taken from the review of W. Glockle et al. They 
are expressed via spin 2x2 matrices (Xj for the nucleon and 3x3 matrices Vi and Vik for 
the deuteron. The latter are related to the deuteron spin matrices Sf 



V2 



1 

1 1 
1 



S-ii 



1 

V2 



-i 

1 -i 
i 



, Sz 



1 

0-1 



(33) 



One has Pi = Si, Vik = 3/2{SiSk + SkSi), Vzz = 3SzSz-2I, andVxx-Vyy = 3{SxSx- SySy). 
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Nucleon analyzing powers Ak are 

A, = IL^. (34) 

ir [aa') 

If the scattering plane is the xy plane and the y axis points to the direction q x q' then due 
to parity conservation = Az = and the only non-zero component is Ay. 
The deuteron vector and tensor analyzing powers are defined as 

_ Tr (aP,at) _ Tr (aP.fcgt) 

~ Tr(aat) ' ^^'^ " Tr (aat) ' ^^^^ 

Parity conservation puts A^jA^, A^y and Aj^^ to zero. So the non- vanishing and independent 
analyzing powers are defined by 

"x/s 1 11 

Also spin transfer coefficients are given in the review. They have the same structure as 
the quantities above, with slightly different matrices to be inserted between a and a) . 



III. RESULTS 

Our results for the differential cross section and nucleon analyzing power Ay (Fig. [T]), 
deuteron vector iTn and tensor analyzing T20 powers (Fig. [2]), and T21 and T22 (Fig. [3]) for 
nd elastic scattering at 3 MeV using the AV14 NN potential are shown together with the 



benchmark calculations of Kievsky et al. 24 1. 



The theoretical predictions are compared with the experimental nd Ay data at 3 MeV 



251]. In both calculations all values of the total three-body angular momentum up to M 



15/2 have been used. In our calculations the total angular momentum of the pair of nucleons 



j23 has been taken up to 3, while in 2J] this value was taken up to j23 = 4. It should be 
noted that in the case of nd scattering increasing 723 by unity raises the number of partial 
waves from 62 up to 98. This difference presumably explains minor differences between these 
two calculations around the maximum values of Ay (Fig. [2]) and of iTu (Fig. [3]), where the 



predictions of Kievsky et al. [241 are consistently higher by about 2-3 %. Differences in T20, 
T21 and T22 are even smaller, about 1 %. 

For the pd elastic scattering at 3 MeV, results of our calculations for the differential cross 
section and proton analyzing power Ay are shown in Fig. H] together with those from the 
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benchmark calculations of Deltuva et al. 26| . Our calculations have been performed using 
the AV14 NN potential and involving the correct asymptotic condition to take into account 
the Coulomb interaction while those of Ref. 
screening and normalization procedure for the Coul 
are compared with the experimental data of Ref. 
angular momentum up to M = 15/2 are used in our calculation, while in Ref. |26| value of M 
is much larger. We chose values of J23 up to 4 (up to 152 partial waves taken into account), 



26i used the AVI 8 NN potential and the 
omb force. All theoretical calculations 



27l |. All values of the total three-body 



whereas in Ref. 



26| these values up to 5 have been used for the strong interaction (207 



partial waves were taken into account). Again this truncation results in a smal 



between our predictions for polarization observables and those from Ref. 26|]. The results 



disagreement 



of calculations for the deuteron vector zTn and tensor T20 analyzing powers as well as the 
experimental data [27| are shown in Fig. [51 The results of calculations for the deuteron 
tensor T21 and T22 analyzing powers as well as the experimental data 27|] are shown in Fig. 
[6l Predictions of our calculations and those of Ref. 26|] are in reasonable agreement. 

In addition to our new results for nd and pd elastic scattering we would like to present 
our new results for pd breakup scattering at Eijjb=14.1 MeV obtained with the Malfliet- 
Tjon (MT) I-III potential. In our paper Ref. |22| results for inelasticities and phase shifts 
were obtained in s-wave approximation for both the strong and the coulomb interactions. 
This means that only partial waves with / = were taken into account for nuclear and 
electromagnetic forces. It was explicitly pointed out and clearly explained in the paper. In 
Table [Jour old and new results together with those of Ref. 26| are given. Our new results 
presented in rows 2 — 5 were obtained for s-wave (MT) I-III potential but the Coulomb 
interaction was taken into account for different choices of sets of basis states (number of 
FNNM equations) in dependence of the maximum value of the two-body angular momentum 
j23- For the calculation we have neglected by the contribution of the basis states with the 
total isospin T=3/2. This negligence makes the relative error less than 0.2%. In the Table 
[J one can see convergence of our results to those from Ref. 26|]. Disagreement between 
inelasticity parameters is about 1% and is about 0.1 degree for phase shifts. As is pointed 



out in Ref. 



26|, the authors used the perturbation method. Our calculations have been 



performed by direct solution of the FNNM equations reduced to a set of linear equations 
with the resulting matrix having tri-block-diagonal structure. Small disagreements between 
results for s-wave pd breakup scattering one can explain by different numbers of partial 
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waves taken into account (up to 126 in our calculations and up to 398 in calculations of 
Deltuva et al.). The authors in Ref. [26| have emphasized that for such large set of basis 
states direct solution is impossible and one has to apply the perturbation theory. 



IV. DISCUSSION 



Our results for nd elastic scattering at 3 MeV and those from the KVP and momentum- 
space calculations are in very good agreement and minor differences can be related to smaller 
values of j23 taken into account in our calculation. In the energy region from 1.2 to 10 MeV 
28| theoretical predictions are 25-30% lower than the experimental data. 

For pd elastic scattering excellent agreement within 1% between momentum-space and 
coordinate-space calculations based on a variational solution using a correlated hyperspher- 
ical expansion predictions at 3, 10 and 65 MeV incident nucleon energies have been demon- 
strated in Ref. j26|. Predictions of our calculation and that of Deltuva et al. 26 1 differ in 



the use of the NN potential and in values of the total three-body angular momenta M and 
of the total angular momenta of the pair of nucleons j23 taken into account. Our prediction 



is about 5% lower for Av than that of Ref. 



26|, and surprisingly about 10% higher for 



iTii. Predictions for tensor analyzing powers agree to better than 5%. Comparison with 
the experimental data of Ref. 27| confirms the Ay and iTu puzzles. Both our calculations 



and those of Ref. 



26| are lower than measured values of Ref. 27(]. Results using the AVIS 



NN potential give better agreement with experimental data for T20 and T21. However, sur- 
prisingly our calculation using the AV14 is in better agreement with the analyzing power 
T22, possibly indicating differences between AV14 and AVIS potentials. 

To end the discussion, we would like to compare our results for polarization observables 
with those from Ref. 18|. In that paper the authors have performed a detail comparative 



study of modern three-nucleon models together in conjunction with the AV18 NN potential 
to calculate observables for pd elastic scattering at Eiab=3 MeV. The authors have shown 
that only the N^LOL TNF model allows to improve the description of Ay and iTu noticeably. 
At the same time the description of T21 becomes slightly worse and there is no change in 
T22- In this regard we would like to note that our predictions obtained with the AV14 NN 
potential and without three-body forces coincide with the experimental T22 data 27| and 



are in good agreement with the result for iTu from Ref. 



18| obtained with three-body 
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forces. 



V. CONCLUSION 

Very good agreement between predictions of our calculations and those of benchmark cal- 
culations demonstrates the soundness of our novel method providing thereby a new approach 
for calculating three-nucleon scattering including nucleon-nucleon and electromagnetic inter- 
actions. Our approach can and will be used to include three nucleon forces and to perform 



additional studies using Kukulin's potential 



and LS modified three-nucleon forces of 



Kievsky [^50,], particularly to study the Ay puzzle. It is well-known that Nd polarization 
observables are the magnifying glass for studying ^Pj forces and calculations that rigorously 
include nuclear and electromagnetic interactions are very valuable. 

Notwithstanding the significance of 3NF, our primary goal is to extend our study using 
AV14 NN potential and including the Coulomb potential to energies above the two-body 
threshold and to focus on breakup data and on established discrepancies. Our next step 
is to use the AV18 NN potential. As discussed in this article, we have already established 
interesting differences in T20, T21 and T22 most likely due to difference between AV14 and 
AV18 NN potentials. 
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TABLE I. pd quartet inelasticity and phase shift (in deg.) at £';aB=14.1 MeV. Nst is the total 
number of coulomb partial waves in the case T=l/2^ For strong interaction only one partial wave 
(Z=0) is taken into account. The results from Ref. [22] obtained with a single coulomb partial wave 
{I = 0) are given in the first row. New results obtained for various sets of basis states are given in 



rows from 2 to 5. 


The result obtained in Ref. 


26 1 using 


perturbation method is given in the last 


row. 
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FIG. 1. Differential cross section and neutron analyzing power Ay for nd elastic scattering at 
3 MeV lab energy as function of the cm. scattering angle. The solid lines correspond to our 
results obtained with AV14 NN potential. The dashed lines correspond to results of Kievsky at al. 
obtained with AV14 NN potential 2j]. The experimental data are from Ref. [25i] . 
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FIG. 2. Deuteron vector iTn and tensor T20 analyzing power for nd elastic scattering at 3 MeV 
lab energy as function of the cm. scattering angle. The notations are the same as in Fig. [TJ 
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FIG. 3. Deuteron tensor analyzing powers T21 and T22 for nd elastic scattering at 3 MeV lab 
energy as function of the cm. scattering angle. The notations are the same as in Fig. [TJ 
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FIG. 4. Differential cross section and proton analyzing power Ay for pd elastic scattering at 
3 MeV lab energy as function of the cm. scattering angle. The solid lines correspond to our 
results obtained with AV14 NN potential. The dashed lines correspond to Deltuva at al. results 
26l | obtained with AV18 NN potential. The experimental data are from Ref. [27l] . 
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FIG. 5. Deuteron vector iTn and tensor T20 analyzing power for pd elastic scattering at 3 MeV 
lab energy as function of the cm. scattering angle. The notations are the same as in Fig. HI 
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FIG. 6. Deuteron tensor analyzing powers T21 and T22 for pd elastic scattering at 3 MeV lab 
energy as function of the cm. scattering angle. The notations are the same as in Fig. HI 
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